Multiple membrane cavity optomechanics 
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We investigate theoretically the extension of cavity optomechanics to multiple membrane systems. 
We describe such a system in terms of the coupling of the collective normal modes of the membrane 
array to the light fields. We show these modes can be optically addressed individually and be cooled, 
trapped and characterized, e.g. via quantum nondemolition measurements. Analogies between this 
system and a linear chain of trapped ions or dipolar molecules imply the possibility of related 
applications in the quantum regime. 
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Cavity optomechanics is an emerging field at the 
boundary between quantum optics and nanoscience. Re- 
sulting in part from experimental innovations at the 
mesoscopic scale, optomechanical systems - mechanical 
systems that can be manipulated by light - have re- 
cently generated much experimental and theoretical in- 
terest P, i, i, 0, @, Q. They offer the prospect of real- 
izing quantum effects at a macroscopic scale Q , of sup- 
plying novel quantum sensors for applications ranging 
from single molecule detection 0] to gravitational wave 
interferometry @, for the quantum control of atomic, 
molecular and optical systems (loj . and f or p ossible new 
quantum information processing devices . 

As cavity optomechanics begins to mature as a field, we 
recognize that scalability is an important aspect of any 
technology. In particular, it is immediately relevant to 
possible uses in information processing. This has been ac- 
knowledged in proposals for constructing quantum com- 
puters using trapped ions, cavity quantum electrodynam- 
ics, neutral-atomic lattices, nuclear magnetic resonance, 
spintronics, dipolar molecules, etc., see Ref. [12] and ref- 
erences therein. It is therefore important to investigate 
the scaling of current cavity-based techniques to a larger 
number of optomechanical elements. 

In this Letter we consider the optomechanical cooling 
and trapping of a small array of partially transparent di- 
electric membranes inside a high-finesse cavity driven by 
laser radiation (Fig. [T]) . Presenting first the case of two 
moving membranes, N = 2, where explicit analytical re- 
sults are readily obtainable, we show that it is described 
most conveniently in terms of the optical cooling, trap- 
ping and measurement of the normal modes of the linear 
chain formed by the membranes. This demonstration de- 
pends crucially on the fact that the optical spectrum of 
the cavity is symmetric in the symmetric normal modes of 
the membrane array, while it docs not have well-defined 
symmetries with respect to the motion of the individual 
membranes. We then extrapolate these considerations 
to the case of N membranes in a cavity, and draw an 
analogy to chains of trapped ions (la UM and of dipo- 
lar molecules [Hj]. We conjecture that this analogy may 
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FIG. 1: A high finesse optical cavity with two mirrors fixed 
at ±3L, and two dielectric membranes centered at qi ~ — L 
and Q2 ~ L respectively. Pi„ signifies the power in the laser 
light that may be used to drive the cavity. 



have implications for information processing if, as ex- 
pected, the membranes can be experimentally placed in 
the quantum regime. 

Our immediate motivation derives from a pioneering 
experiment that demonstrated that a single membrane 
can be optomechanically cooled when placed inside a 
high finesse optical resonator This is the first ex- 

periment to show that the two technologically challeng- 
ing requirements of high optical and mechanical quality 
can be allocated separately to the end-mirrors and the 
dielectric membrane respectively, distinct from the stan- 
dard two-mirror cavity configuration This work also 
pointed out that in addition to the usual linear coupling 
of radiation to the displacement q of the membrane that 
facilitates back-action cooling, a coupling proportional 
to q 2 is also possible. If the membrane oscillation fre- 
quency is larger than the optical linewidth of the cavity 
[l7 . la . 19j |. this quadratic interaction enables quantum 



non-demolition measurements of the energy of the vi- 
brating membrane. Linear couplings do not allow such 
measurements [20| . 

The present work draws on our previous analysis of 
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the case of a single membrane [2l|, which yields an ap- 
proximate microscopic Hamiltonian that describes the 
optomechanics of the system. We show that this for- 
malism can in principle be generalized to TV membranes, 
with the generic features extracted simply from consid- 
ering the case N = 2 (Fig.[T|). More detailed calculations 
will be presented elsewhere [22[ . 

Our starting point is a cavity with two fixed and per- 
fectly reflecting end mirrors and two identical vibrating 
dielectric membranes, each of reflectivity R, mass m and 
mechanical frequency uj m , see Fig.[T] We model the sys- 
tem with the quantum mechanical Hamiltonian 21 1 
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fiwi(qj)a\au (1) 



where pj and qj, j = 1, 2 denote the momentum and po- 
sition of the membranes and Wj, i ~ 1,2,3, are the reso- 
nance frequencies of three near-degenerate optical modes 
with creation and destruction operators and a\ respec- 
tively. The commutation rules obeyed by these operators 



are [<7j,pz] = ihdji and [c 



We discuss shortly 



why three optical modes are sufficient for a minimal de- 
scription of the system. 

The first term in Eq. (QJ is the mechanical energy of 
the oscillating membranes, and the second term the en- 
ergy of optical modes of the full resonator of length 6L. 
The dependence of the cavity mode frequencies w.j on the 
positions qj of the membranes is central to the descrip- 
tion of the system, since it determines the optomechan- 
ical couplings [2lj]. The u>i are obtained by solving the 
classical Maxwell equations for the full resonator. We as- 
sume that the moving membranes are much thinner than 
an optical wavelength and model them by spatial delta 
functions [2l[. For R = 1, the resonator consists simply 
of three uncoupled cavities whose eigenfrequencies in the 
absence of membrane motion (yi.2 = TL) are threefold 
degenerate and are given by 



mrc 
~2L' 



(2) 



Here n is a positive integer, c is the velocity of light, and 
2L is the length of each sub-cavity. This is the reason we 
included only three modes in Eq. |T]); the spectrum may 
generally be grouped in such triplets. 

For R ^ 1, the three resonators are coupled, this cou- 
pling lifting the degeneracy of the frequencies u> n . Solv- 
ing Maxwell's equations with the appropriate boundary 
conditions at the membranes, and neglecting light ab- 
sorption in these membranes, we find that the allowed 
wave-numbers of the complete resonator can be expressed 
in terms of the relative coordinate q = q\ — q2 and the 
center-of-mass (COM) coordinate Q = (qi + q2)/2 of the 
membranes via the trigonometric equation 



sin 2(6 + 3kL) + sin 2 6 sin 2k (3L - q) 
= 2 sin 9 cos (8 + kq) cos 2kQ, 



0>1,3 
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FIG. 2: A portion of the optical frequency spectrum for the 
two-membrane cavity. A closely spaced triplet of frequencies 
= 1,2,3) is shown, using a higher resolution near the 
mode equilibrium point go = 2. For Q / the optical mode 
frequencies are similarly modulated in a direction orthogonal 
to q. The cavity parameters used in the plot are R — 0.5, L = 
1. 



where sin 9 = yR. 

Equation ([3]) is invariant under the transformation 
Q — ► —Q. Physically this means that the frequency 
spectrum is symmetric in the COM motion. We show 
below that this symmetry can be exploited to obtain a 
description of the system purely in terms of the 'phonon' 
modes associated with the relative and COM coordinates 
q and Q separately, that is, excluding any effects result- 
ing from the coupling between these coordinates. This is 
in contrast to an equivalent description in terms of the in- 
dividual oscillator coordinates q\ and q2, which typically 
retains coupling terms. 

Equation ([3]) can be solved numerically for a given set 
of system parameters. Figure [2] shows a triplet of optical 
modes u> n ^(q, Q) in the vicinity of the equilibrium point 
of relative motion qo = 2 and for Q = 0. The optical 
frequencies are strongly modulated along q, a situation 
familiar from the case of a single membrane, N = 1. 
These modulations arise from the avoided crossings asso- 
ciated with the lifting of the degeneracy of the frequency 
triplet due to the coupling of the sub-cavities [13] ■ Simi- 
lar periodic modulations appear in u> n ,i(q,Q) for fixed q 
and Q varying. 

Considerable insight can be gained from approximate 
analytical solutions of Eq. (J3j) [22| . Perturbation theory 
in the small parameters (q — qo, Q — Qo) <C 2irc/u! n yields 
for example the triplet (i = 1,2,3) of optical frequencies 
in the vicinity of any go and Qo 

w n ,i(«,Q) = &n,i + B» ti (q - qo) + B'i(Q - Qo) 



+M nii (q - qof + M' n<i (Q 

+P n ,i{<i - qo)(Q - Qo) + ■ 



Qo) 2 (4) 



(3) 



The various terms in this equation fully describe the 
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basic optomechanical properties of the system. Specif- 
ically, A n ^ are the frequency shifts due to the sub- 
resonators coupling in the absence of membrane motion. 
B n i and B' n i determine the strength of the linear op- 
tomechanical couplings associated with the "breathing" 
and COM modes of motion, respectively, producing a 
back-action of the mirror motion on the light field that 
can be exploited in mirror cooling. M n> i and M[ gov- 
ern the quadratic couplings and can lead to quantum 
non-demolition (QND) energy measurements |l6[). Fi- 
nally P n> i quantifies the coupling between the relative 
and COM modes and is responsible for normal mode de- 
coherence as well as down-conversion [23|. We discuss 
later in this Letter circumstances under which that cou- 
pling can be cancelled. 

Consider for concreteness the case qo = 2L, Q = 0. In 
this case we have 
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sin 



and 



f n ,i sin(6'„ !l )sin(26') 



3 cos 



i(0„,i) sin 2 (6>) + 3 cos(26> + 30^) + sin(20) sin(0 n>i ) ' 



(6) 



where 6 n>i = 2A lM L/c and £, M = 6 n ^/(2L) 2 . 

In that case, both B n ^ and B n ^, which are the slopes 
of the two upper curves in Fig.[2j are different from zero. 
The linear spatial dependence of the optical frequencies 
on q can therefore be used for back-action cooling and 
trapping of the "breathing" mode of the pair of mem- 
branes. In contrast, -B ra> i = 0, and the corresponding 
resonator mode cannot be used for that purpose. How- 
ever it can be used for a non-demolition measurement of 
the q mode, as it provides a purely quadratic coupling 



M n ,x =--T^,itan0^O, 



(7) 



where r = AL/c is the round trip time for each sub- 
cavity. Hence both linear and quadratic optomechanical 
couplings are accessible at the same point qo = 2 using 
modes from the same frequency triplet. This is unlike 
the situation for a single membrane, see Refs. [3, EH. 

The analysis for the mode Q is similar to the case of a 
single membrane since both describe the effects of COM 
motion [2l[. The choice of Qo = implies the absence 
of the terms linear in Q in the expansion of Eq. ((3]), a 
direct consequence of the fact that Eq. ([3]) is symmetric 



about Qo = 0. Hence B' ni — 0. The quadratic cou- 
pling is however non-zero; for reasons of space we will 
provide an expression for it elsewhere |22j . On the other 



hand, when the membranes are positioned so that Qo ^ 
then non-vanishing terms linear in Q arise in Eq. ^ and 
optomechanical center-of-mass cooling becomes possible. 
Further analysis of this case will be presented elsewhere 
0. 

We finally turn to the last term in Eq. ((4|). As al- 
ready mentioned, it describes the coupling between the 
breathing and COM modes of motion of the membranes. 
In perturbation theory this bilinear term is smaller than 
the linear terms, hence it does not effect significantly the 
back-action properties used in cooling. Furthermore, for 
Qo = we must have that P n ,i = since Eq. ((3]) is even 
in Q and hence the expansion cannot contain terms that 
are odd in Q. To lowest order, the coupling between these 
modes of motion is then proportional to (q — q )Q 2 and is 
negligible. In that case, the breathing and COM modes 
remain to an excellent approximation the 'true' normal 
modes of the full two-membrane optomechanical prob- 
lem. This is a central result of this paper. In physical 
terms, this implies that with an appropriate combination 
of optical frequencies, it is possible to independently cool, 
trap and characterize the breathing and COM modes of 
the 2-membrane array. 

We now discuss in broad terms the extension of these 
considerations to an array of iV > 2 membranes. The 
addition of a membrane positioned at a node of the in- 
tracavity field does not change the optical finesse mea- 
surably ; if the mirror is placed away from a frequ ency 
extremum, the finesse is multiplied by R ~ 0.99 3, 2lf - 
Thus a scaling up to N ~ 10 membranes should be prac- 
tical. 

The boundary conditions of the optical resonator yield 
a trigonometric equation similar to Eq. (|3|) independently 
of the number N of membranes. It is always possi- 
ble to express this equation in terms of the N normal 
modes Qj of the array, since they are linear combina- 
tions of the individual oscillator coordinates qj. Gen- 
erally the resulting optical spectrum will be periodic 
in every mechanical mode, guaranteeing the existence 
of ranges where the dependence of the N + 1 frequen- 
cies uj nt i({Q j}) , i = 1, ...N + 1, exhibits extrcma, and 
other regions of linear dependence. This in turn implies 
the possibility of cooling, trapping and performing QND 
measurements of the energy of all collective modes. The 
optomechanical parameters coupling any optical mode 
to any mechanical mode can be extracted either analyti- 
cally, from expressions such as Eqs. (0][7|), or numerically. 
The cooling and trapping effects due to back-action , be- 
ing linear in the Qj, are always independent of the cou- 
pling between the phonon modes. 

Most importantly, we expect the optical spectrum to 
be symmetric in every symmetric normal mode of the ar- 
ray. We have shown this explicitly for the COM mode in 



4 



the cases N = 1 [2lJ and N = 2 (this work). There 
are N/2 symmetric modes for N even, and N/2 + 1 



for N odd [24|. For N = 3 for example, with indi- 
vidual oscillator coordinates qj,(j = 1,2,3) we expect 
the optical spectrum to be symmetric in the COM mode 
Qi = (<Zi + 92 + <?3)/3 as well as in the 'scissors' mode 
Q2 = (<7i — 2q2 + q,s)/6, but not in the stretch mode 
Q3 = (<73 — [22|. However the amount of symmetry 

available is adequate to cancel all mode-mode coupling 
terms by appropriate positioning of the membranes. For 
N > 4, some mode coupling terms will have to be in- 
cluded in the description of the optomechanics. For N 
even, for example, ( N l 2 ) = N(N — 2)/8 coupling terms 
out of a total possible (^) = N(N — l)/2 will remain in 
the expression for each frequency. 

The full behavior of the system may be modelled by 
a Hamiltonian analogous to Eq. (Q]), combined with a 
quantum treatment of the noise associated with the in- 
put fields and with dissipation, see e.g. Ref. [l8j|. Exper- 
imentally the noise spectrum of any collective modes can 
be obtained in principle by modulating the light field at 
the specific mode frequency 0. 

Our description of the membrane array is reminiscent 
of a chain of trapped ions [l3|, [l4 , 23, 24j, with however 
significant differences. First, in the present case the cou- 
pling between mirrors is due to radiation pressure rather 
than the Coulomb interaction, and is therefore switch- 
able. As such, the situation is perhaps more akin to the 
case of dipolar molecules, where the coupling is due to 
the dipole-dipolc interaction and can also be tuned. Also, 
membranes and mirrors do not possess internal degrees 
of freedom, at least at the level of the present descrip- 
tion. Still, the analogy to the physics of trapped ions 
or dipolar molecules on lattices raises the possibility of 
using multi-membrane systems as information processing 
devices. Each membrane can be addressed individually 
through its suspension, and additional degrees of freedom 
such as rotation could possibly be added and exploited 
f25l | . Issues of noise and decoherence will of course be ex- 
tremely challenging, but some work along these lines is 
already being carried out in efforts to cool such systems 
to their vibrational ground state 17, [lij. In addition, 
as is the case for trapped ions the independent control 
of the collective normal modes of the array will become 
technically increasingly challenging with increasing N. 

In conclusion we have considered the extension of 
cavity optomechanics to multi-membrane systems, and 
found that the normal modes of the system can cooled, 
trapped and measured. The similarity of this system 



to chains of trapped ions or dipoles leads us to conjec- 
ture that it may find possible future application as an 
information-processing device operated in the quantum 
regime. 
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